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Rigorous bounds are established for the expansion coefficients governing the shape of semileptonic 
K — > 7r form factors. The constraints enforced by experimental data from r — > K-kv eliminate 
uncertainties associated with model parameterizations in the determination of |V U s|. The results 
support the validity of a powerful expansion that can be applied to other semileptonic transitions. 
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INTRODUCTION 



Semileptonic K — » -kIv (Ka) decays provide the 
most precise determination of the Cabibbo-Kobayashi- 
Maskawa (CKM) matrix element \V US \ Presently, the 
dominant uncertainty in the experimentally determined 
quantity | V^ s ^+ (0) | arises from the unknown shape of 
the hadronic form factor F+Ja 2 ), as a function of mo- 
mentum transfer q @, H, 0, Q. The data can be fit to 
both simplified pole models and series expansions, but 
the uncertainty inherent in these simplifications is diffi- 
cult to estimate, and the difference between the resulting 
|y us F + (0)| determinations represents a systematic error. 
Constraining this shape from first principles and provid- 
ing rigorous error estimates is an important problem. 

The remaining (and presently, dominant) error in the 
K a determination of | V us | arises from the normalization, 
F + (Q). K — > ntv data (specifically K^) also constrain 
the shape of the scalar form factor, Fo(q 2 ). This can po- 
tentially reduce dominant errors in the theory normaliza- 
tion, either by avoiding an extrapolation from zero recoil 
in lattice determinations Q , or by fixing low-energy con- 
stants of chiral perturbation theory [tJ. At present, the 
comparison between theory and experiment, and between 
different experiments, is complicated by uncertainties in 
the form factor parameterization. Again, extracting as 
much information as possible from the experimental data 
without model assumptions is an important task. 

This paper establishes rigorous bounds for the expan- 
sion coefficients appearing in a general parameterization 
of the semileptonic K — > 7r form factors. The framework 
for the analysis is based on familiar arguments invoking 
analyticity and crossing symmetry Q. However, in con- 
trast to conventional arguments appealing to unitarity 
and the evaluation of an operator product expansion 
(OPE), bounds for the vector form factor are derived 
directly from experimental data for hadronic tau decays, 
r — > K~kv. The result is more stringent than can be 
obtained from the OPE analysis, and applies to a more 
general class of parameterizations. 



Analyticity and convergence. Form factors are defined 
as usual by the matrix element of the weak vector current 
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F+{t — q 2 ) and Fo(t) can be extended to analytic func- 
tions throughout the complex t plane, except along a 
branch cut on the positive real axis starting at Kir pro- 
duction threshold t = t+ [t± = (mx ± m^) 2 ]. The cut 
plane is mapped onto the unit circle by 
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where to € (— oo,£ + ) is the point mapping onto z — 
0. [22| The form factors are analytic in \z\ < 1, and so 
may be expanded in a convergent power series: 
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where (f> is an as- yet arbitrary function analytic in | z \ < 1 , 
which may depend on one or more parameters, denoted 
generically (and for reasons to become clear) by Q 2 . 

The function z(t, to) sums an infinite number of terms, 
transforming the original series, naively an expansion in- 
volving t/t + < 0.3, into a series with a much smaller 
expansion parameter. For example, the choice to — 

(1 — y/l — t-/t+) minimizes the maximum value of z 
occurring in the semileptonic region, and for this choice 
\ z {t, to) | 0.047. The function cj> and the number to m ay 
be regarded as defining a "scheme" for the expansion. 
The expansion parameter z and coefficients afc are then 
"scheme-dependent" quantities, with the scheme depen- 
dence dropping out in physical observables such as F(t). 

Neglecting terms beyond z N in (|3|) introduces a rela- 
tive normalization error AF/F = 0(z N+1 ). Similarly, 
the error on the relative slope involves terms of order 
(N + l)z N . Since \z\ 2 < 2 x 10~ 3 and |z| 3 < 1 x 10~ 4 , 
simple power counting in z yields strong constraints 
on the impact of higher-order terms in the expansion, 
provided that the coefficients ak/ao are well-behaved. 
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Crossing symmetry and form factor bounds. It is impor- 
tant in practice to determine whether large "order unity" 
coefficients au /a could upset the formal power counting 
in z. To address this question, a norm may be defined 
via 
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By crossing symmetry, the norm can be evaluated using 
form factors for the related process of Kit production. 
The following sections investigate bounds on the integral 
appearing on the right hand side of (|4|). 



II. VECTOR FORM FACTOR CONSTRAINTS 

To compare with unitarity predictions, and to motivate 
a choice of in ([3|) , we consider the correlation function, 



Tl^(q) = i I dW^OlTjV^F^O)*} |0) 

= (-.9 M Y + q^lndq 2 ) + q^Tl (q 2 ) ■ (5) 

An unsubtracted dispersion relation can be written for 
the quantity: (Q 2 = -q 2 ) 



Xi(Q 2 ) = 7T 



i d 2 



l f°° 

[g 2 nr] =~ dt 
~ Jo 



timn^t) 



(G) 



2d{q 2 ) 2 Ly 11 nj {t + Q 2 f 

Assuming isospin symmetry, the contribution of all Kit 
states to the (positive) spectral function Imlli(i) is 
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Choosing: [note that \z\ = 1 along the contour in ^ ] 
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, z(t,t )\- 1/2 ( z(t,t_) y 3/4 t + -t 
1 lo—T ) {l~T) (t + ~t y/± ' (8) 

then yields the inequality: 23] [note that ao(to,Q 2 ) = 
<t>{t ,to,Q 2 )F(t )} 
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For Q 3> Aqcd, Xi(Q ) can be reliably calculated usin 
the OPE in ([5]) . Collecting results from the literature 
[l2j], we have at renormalization scale \i = Q: 24 1 
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FIG. 1: Bounds on the expansion coefficients for the vector 
form factor. The top (light) band represents the unitarity 
bound, and the lower (dark) band is a direct evaluation from r 
decay and perturbative QCD. The perturbative contribution 
is shown separately as the dashed line. 
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FIG. 2: Bounds on Af + in as a function of to. The 
dashed lines are for the default choice of 4>f + in (JSJ, with 
Q = (top) and Q = lOGeV (bottom). The solid band gives 
the (Q-independent) result when (/>f + = 1. 



where corrections of order m u /m s are neglected, 
and m s (2GeV) = 0.087(8) GeV O, -m s (uu) - 
a s (G 2 )/12n = -0.0001(8)GeV 4 [3. The light band 
in Fig. [1] shows the resulting bound on the quantity 
A F+ (to,Q 2 ), setting t = and F+(0) w 1. The per- 
turbative uncertainty is estimated by varying fj? from 
Q 2 /2 to 2Q 2 , and allowing for higher-order contributions 
of relative size ±1 X a*. Uncertainties from perturbative 
and power corrections are small above Q — 2 GeV, but 
become significant below this scale. The width of the 
band represents a la contour obtained by adding uncer- 
tainties in perturbative and power corrections linearly. 
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Although the OPE breaks down at small Q, the norm 
in Q remains perfectly well defined. The dark band in 
Fig.Q]shows the result of evaluating the integral in (j4|) us- 
ing t — > Kttv decay data [l5[ for the region t + < t < m 2 ., 
and estimating the remaining contribution from t > m 2 . 
using perturbative QCD. Uncertainties from the r decay 
data are conservatively estimated by taking the maxi- 
mum (minimum) value of the weight function in ^ for 
each bin, and multiplying by the corresponding upper 
(lower) bound of the la interval for the Kir component 
of the spectral function measured in [l5T ] . It may be noted 
that for the present purpose, there is no need to resolve 
the underlying resonance structure of the spectral func- 
tion, nor to subtract the (small) 5-wave contribution. 

The remaining perturbative contribution is estimated 
using F+(t) ~ iirf-wfKOis/t [16(, conservatively setting 
a s w 1, and assigning 100% uncertainty to the result. 
The effect is AA 2 F+ < 2 for all values of Q e (0, 10 GeV), 
and is insignificant when added in quadrature to the 
much larger contribution from the resonance region. 

As illustrated by Fig. [21 the dependence of Ap + on 
t is very mild. Here we again set F + (to) ~ 1; this is 
strictly an overestimate, and accounts for much of the to 
variation in the figure. The norm for all schemes © with 
< t Q < t- and < Q < 10 GeV satisfies 



A F+ (t ,Q 2 )<13. 



(11) 



More refined values for Ap + can be obtained from Figs.[T] 
and Oat particular values of to and Q. 

The result from the OPE represents an upper bound for 
Ap + . For large Q, the leading contributions to Ap + 



are 
n2 



O(Q ), coming from the regions t ~ Aq CD , and t 
in J3J). [H[ The OPE overestimates A F+ by a factor 
~ (Q/Aqcd) 2 - As Fig.[T]shows, although the OPE evalu- 
ation becomes exceedingly precise, the resulting unitarity 
bound begins to wildly overestimate the norm of the form 
factor already at Q ~ 2 — 3 GeV. 



III. SCALAR FORM FACTOR CONSTRAINTS 

In the scalar case, an unsubtracted dispersion relation 
can be written for the quantity 
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FIG. 3: Bounds on the expansion coefficients for the scalar 
form factor. 



yields the inequality 
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Again, Yn(Q 2 ) can be reliably calculated when Q ^> 
Aqcd [SQ: 
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where m s (2 GeV) = 0.087(8) GeV p] 



(ss)) 



-0.0006(3) GeV 4 [lj. The light band in Fig. shows 
the la contour obtained by adding errors from all sources 
linearly, as described after (fit)]) : in this case, higher-order 
contributions of relative size ±10 x a 4 are included. [26[ 
The behavior of the unitarity bound with respect to Q 
is milder than in the vector case. For large Q, the leading 
contributions to Ap are O(Q ), coming from the regions 
t — Aq CD and t — Q 2 in (gj). [27| Again, it should be 
remembered that the OPE gives an upper bound to Ap , 
although for this case the bound is overestimated by only 
a single power of Q /Aqcd ■ While it is not possible to rig- 
orously bound the norm of Fo using the OPE for schemes 
with Q < Aqcd, there is no parametric enhancement of 
the integral in (fTS"|) for these Q values. A conservative 
bound on the norm for all schemes with < to < t— and 
< Q < 10 GeV is 



A Fo (t ,Q 2 ) <5. 



(17) 



In principle, this bound could be improved by extracting 
the Kir component of the scalar spectral function from 
r — > Kirv. 
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IV. IMPLICATIONS FOR K a AND \V US \ 

To implement the expansion J3]) we must choose a 
"scheme" specified by the functional form of (j), and by 
the values of Q 2 and to. Possible candidates for <p are 
= 1, or (j) given by ((8]) and (fT4|) . The latter choice en- 
sures that the norm in (TJJ is independent of to and relates 
in a simple way to the spectral functions extracted from 
r decay; with, say Q 2 = 2GeV 2 , this choice also allows 
comparison to OPE bounds. [28[ 

The value of to can be selected according to various cri- 
teria. For example, the choice to = t+(l — y/1 — t_/t+) 
minimizes the maximum value of z throughout the 
scmilcptonic range, while to = simplifies the transla- 
tion between coefficients a k and the derivatives of the 
form factor at t = 0. Alternatively, the choice of to can 
be used to minimize correlations between shape parame- 
ters. 

When only a single form factor contributes (as in K e ^ 
decays), the error correlations between expansion coeffi- 
cients dk(to,Q 2 ) measured in an experiment with near 
ideal acceptance and resolution are completely deter- 
mined without a priori knowledge of the form factor 
shape. With Ar,/r 4 cx l/y/Tl where Ti ~ p(U) 3 F + (ti) 2 
are partial rates measured in intervals of t, it follows that 
the error matrix U determined from a x 2 fit is: (here 
p = \p\ is the pion momentum in the kaon rest frame) 

Ukk ' - 2 da k da k , * Jo 0f + (t, to, Q 2 ) 2 ^ ' 

(18) 

If the K e 3 data is fit to a normalization proportional 
to F + (0) and two shape parameters a\(to, Q 2 )/ao(to, Q 2 ) 
and 02(^0, Q 2 )/ao(to, Q 2 ), the correlation between ax/do 
and ai j cio vanishes in the ideal case for to w 0.39 1_ (with 
F+ as in © and Q 2 = 2GeV 2 ). [H 



A. F+ and AT e 3 

Only the vector form factor is relevant for the elec- 
tron mode. From Jill) and \z\ < 0.05 [at t = t + (l - 
i/l — t-/t + )], the quadratic term, a,2Z 2 in ([3]), can af- 
fect the normalization by an amount not greater than 
^4f + M 2 < 0.029. To establish a precision significantly 
smaller than this, ai must be constrained by data. The 
cubic term is bounded by Ap + \z\ 3 < 0.0014. These errors 
can also be robustly estimated by including higher order 
terms in fits to the experimental data, under the con- 
straint (fTT]) . Extraction of other observables, such as the 
form factor slope and curvature at t = can be treated 
similarly. 

The simple Taylor expansion in t converges much more 
slowly than the z expansion in ([3]), and it is difficult to 
reliably bound the errors associated with higher-order 
terms in fits to the truncated t series. The systematic 
approach advocated here should remedy this situation, 



and make a more definitive comparison between different 
experiments possible. [3(| 



B. F and K„ 3 

The scalar form factor can be probed in the muon 
mode. With enough precision, constraining the slope of 
this form factor could aid lattice and chiral perturbation 
theory estimates of F + (0) = F o (0) by eliminating extrap- 
olation errors [g], or by determining higher-order con- 
stants in the chiral Lagrangian Q- In a general scheme: 
(the prime denotes derivative with respect to t) 

Fq-(O) ^,,(0, to, Q 2 ) , ^ + ^(0,fd) + ... 
Fo(0) M(0,t ,Q 2 ) lU ' toJ l + £z(0,t«,) + ... ' 

(19) 

Present experiments do not strongly constrain the coef- 
ficients beyond the linear a\z term in ([3]). From (fT?)) 
and \z\ < 0.05 it is straightforward to estimate the re- 
sulting uncertainty on the slope: A[m 2 i ;l Q(0)/Fo(0)] w 
Ap a x 0.001 < 0.005. This error can also be robustly 
estimated by including higher order terms in fits to the 
experimental data, under the constraint (|17[) . 

To further constrain the slope requires either more pre- 
cise semileptonic data, or independent constraints on the 
form factor. One type of constraint appeals to symme- 
try principles, in particular at the Callan-Treiman point 
where Fo(Aj<- ff ) =_£k / fn+Acr with a controllably small 
correction, Act UM- Unfortunately, this point is too far 
removed from t = to have a very strong influence on the 
extraction of Fg(0)/Fo(0). [3l| As a second type of con- 
straint, the phase of the form factor above Ktt threshold 
can be determined from Ktt scattering. Even supposing 
arbitrarily precise scattering data were available up to 
the inelastic threshold, ti ne i = (tuk + 3mj) 2 , it is easy 
to see that dramatic constraints are not expected in the 
semileptonic region. Since the additional data relates to 
the form factor at points along the cut in the t-plane, 
where \z(t, to)| = 1, constraints from the scattering data 
can be satisfied by tuning higher-order coefficients a k in 
([3]), which have relatively minor impact in the semilep- 
tonic region. More quantitatively, with arbitrarily precise 
data, a suitable phase redefinition of the form factor (a 
particular choice of <fi) can postpone the branch point in 
4>F until tinoi- Replacing t + with t- mc \ in , we see that 
the maximum size of the effective expansion parameter 
can then be made as small as \z\ < 0.021 throughout the 
semileptonic region, compared to \z\ < 0.047 when the 
branch point appears at t + . Thus while some improve- 
ment is possible by directly incorporating the experimen- 
tal scattering data, more dramatic gains in precision re- 
quire further dynamical inputs or assumptions [20(. 
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V. IMPLICATIONS BEYOND Kea 

The general parameterization ([3]), with constraints 
(fTTj) and (fT?)) . can be used to systematically analyze the 
K — » 7r£i/ data without model assumptions. This should 
eliminate a dominant uncertainty in the experimental de- 
termination of |V^ s .F+(0)|, and allow for a more definitive 
comparison of shape observables measured by different 
experiments. 

Ki3 decays provide a unique opportunity to investigate 
the convergence properties of the expansion ([3]). Pre- 
cision data can be used to directly constrain the first 
few coefficients, and the existence of a heavy lepton, 
m T 3> rriK + iri-n , makes it possible to establish bounds 



on the dominant vector form factor that are both more 
stringent and more generally applicable than those ob- 
tained from an OPE analysis. This provides a direct test 
of scaling arguments that apply also in the charm and 
bottom systems [10[ . 
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